Highly tunable low-threshold optical parametric oscillation 
in radially poled whispering gallery resonators 
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Whispering gallery resonators (WGR's), based on total internal reflection, possess high quality factors in a 
broad spectral range. Thus, nonlinear optical processes in such cavities are ideally suited for the generation of 
broadband or tunable electromagnetic radiation. Experimentally and theoretically, we investigate the tunability 
of optical parametric oscillation in a radially structured WGR made of lithium niobate. With a 1.04-pim pump 
wave, the signal and idler waves are tuned from 1.78 to 2.5 pirn - including the point of degeneracy - by varying 
the temperature between 20 and 62 °C. A weak off-centering of the radial domain structure extends considerably 
the tuning capabilities. The oscillation threshold lies in the mW-power range. 

PACS numbers: 42.60.Da, 42.65.Yj, 42.79.Nv 



Optical micro-resonators (micro-cavities) attract increasing 
research interest owing to their outstanding properties and 
promises for numerous applications lUrfl. Ultra-high qual- 
ity factors of optical modes, reaching Q s» 10 11 |5[], together 
with small mode volumes provide unprecedented conditions 
for shaping and enhancement of light-matter interactions. Al- 
ready existing and highly demanded applications of these fas- 
cinating new optical elements span from quantum electrody- 
namics and optics to optical filters and sensors HI 51 lit. 

Bringing nonlinear-optical effects to the range of low- 
power continuous-wave (cw) light sources is one of the 
biggest challenges. Whispering gallery resonators (WGR's), 
made of crystalline nonlinear media, are best suited for this 
purpose The quality factors, which are restricted 

from above by light absorption, are among the highest ones. 
Together with a /im-size transversal confinement, a ~ 10 7 in- 
tensity enhancement can be reached. The discrete WGR mode 
structure is well-understood nowadays and the techniques for 
its engineering and coupling light in and out are well devel- 
oped | 
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For most of the nonlinear-optical phenomena, phase match- 
ing is a key issue. In WGR's it acquires specific fea- 
tures: While the most excitable equatorial and near-equatorial 
modes can be treated as plane waves, one has to take into ac- 
count a discrete geometry-dependent mode structure as well 
as a modified dispersion relation [9]. The main impact of the 
discreteness is here in the reduced possibility to meet phase 
matching. The smaller the WGR size, the stronger is this im- 
pact B12I1 . Usually, the fine adjustment to narrow -band cavity 
nonlinear resonances can be made by varying the temperature 
or applying an electrical field 11311 . 

A number of important nonlinear effects have been realized 
in WGR's with low-power cw light sources during the last 
years lEl3l [l3lj20ll . In optical materials, Raman scatter- 
ing and lasing, third-harmonic generation, low-threshold op- 
tical oscillation owing to resonant four-wave mixing, and op- 
tical comb generation via Kerr nonlinearity were reported. In 
X^ 2 ' materials, second-harmonic generation, third harmonic 



generation by cascading two second-order processes, and 
parametric down conversion were demonstrated recently. 

Here we tackle the fundamental problem of tunability of 
nonlinear processes in WGR's by indicating how to vary op- 
tical frequencies over wide ranges. Specifically, we demon- 
strate for the first time a highly tunable low-threshold optical 
parametric oscillation (OPO) in a radially poled WGR made 
of lithium niobate (LiNbOs). 

As it is known from the studies of x^ 2 ) nonlinear-optical 
effects in bulk crystals, quasi-phase matching via domain en- 
gineering filtl provides great scope for tuning J22ll . In the 
WGR case, a radial structure with an even number 2N 3> 1 
of domains, see Fig. la, is best suited for that purpose. For the 
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FIG. 1. Geometric properties of radially poled WGR's: Ideal and 
distorted radial domain structures - a) and b) - and representative 
spectra of the reduction factor \Sj \ - c), d), and e) - for 2A^ = 686 
and SR/R = 0, 0.003, and 0.03, respectively, as defined by Eq. (1). 
Light and dark gray regions in a) and b) indicate the alternating do- 
mains, x is the circumference coordinate, R and r are the large and 
small WGR radii. 
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desired nonlinear process, the optimal large radius R at a spe- 
cific domain number can be evaluated taking into account the 
geometry-dependent corrections to refractive indices of the in- 
teracting waves 1 10]. 

In polar optical materials, including lithium niobate, do- 
main inversion results in changing the sign of the nonlin- 
ear susceptibility x • To account for the 2ttR periodicity 
of the WGR case, we represent the sign-changing distribu- 
tion x (x), where x is the circumference coordinate, by the 
Fourier series 

OO 

S(x)= ]T 5 j G xp(i^), S J= S*, (1) 

j=-oo 

where the Fourier coefficients are given by the average over 
the circumference, Sj = (S(x) exp(— ijx/R)), and S(x) — 
sign[x^ 2 \x)]. The l/R separation between the spatial fre- 
quencies, dictated by the 2ttR periodicity, is indeed the same 
as the wavevector difference for the equatorial modes. As 
soon as the domain pattern (the points of the sign change) is 
known, one can calculate Sj numerically. The quantity | Sj \ 
represents the reduction factor for the nonlinear coefficient 
when using the j-harmonic; it depends solely on the WGR ge- 
ometry. If the domain pattern is strictly periodic, the Fourier 
components Sj can easily be calculated analytically. Only the 
harmonics with j = ±N, ±3iV, . . . are nonzero in this case, 
and the reduction factor for the main harmonic, \Sn\ = 2/tt, 
is close to 1. Any distortion of the periodicity of the domain 
structure makes nonzero all spatial harmonics Sj. 

In our case, the most significant distortion originates from 
off-centering of the radial structure, see Fig. lb. It gives close 
side harmonics with j = iV± 1, NzL 2, etc. Subfigures lc) to 
le) illustrate the impact of off-centering on the Fourier spec- 
trum of the reduction factor \Sj\ for 2N = 686, which is 
relevant to our experiment. One sees that the first side har- 
monics become comparable with the main one already for 
SR/R = 0.003. For SR/R = 0.03, we have a rich dis- 
crete spectrum, where the main harmonic is not the dominat- 
ing one. The reduction factor remains pretty large for many 
spatial harmonics. Failure in the inversion of part of the do- 
mains is not critical for the spectrum, because the long-range 
order is maintained. 

Each Sj-peak can generally be used for phase matching. In 
the OPO case, the corresponding phase-matching conditions 
read 

1 1 _j_ 1 n p n s _|_ n\ _|_ j ^ 

Ap A s Ai Ap A s Aj 2itR 

where A p . s ,i are the pump, signal, and idler wavelengths. 
The effective refractive index n, entering these relations, is 
a known function of A, T, R, and r fldiUM . 

Figure 2 shows tuning curves calculated numerically from 
Eqs. d2j for a pump wavelength A p = 1.04 /im, the WGR 
radii R = 1.54 mm and r = 0.6 mm, extraordinary polariza- 
tion for all three waves, and several values of j. Remarkably, 
there is a substantial difference in the tuning curves even for 



neighboring -peaks. In essence, the tuning possibilities can 
be considerably extended in the non-ideal case when using 
several strongest spectral peaks instead of a single one. This 
extension occurs at the expense of a modest lowering of the 
nonlinear coupling strength leading to a modest increase in 
the threshold intensity. 

According to Fig. 2, phase matching is possible only for 
j > 342 above room temperature. Closer examination shows 
that tuning by heating of the WGR becomes possible only 
owing to the impact of the off-centering if the large radius 
R exceeds noticeably the optimum value of 1.54 mm. Note 
lastly that the points of degeneracy, where A S (T) = Aj(T), 
are almost equidistant for the neighboring tuning curves with 
a temperature step AT w 22°C. 

In order to realize the extended OPO tunability, we have 
fabricated a radially poled WGR with 2N = 686, R ~ 
1.55 mm, and r ~ 0.6 mm. The corresponding R/r ra- 
tio is known to be favorable for coupling light in and out 
with a rutile (Ti02) prism [11]. The resonator was fabricated 
from a 500 fim thick z-cut wafer of congruent lithium niobate. 
The domain structuring was performed using the standard 
electric-field-poling technique with liquid electrodes 12411 : a 
radially patterned photo-resist corresponded to the desired do- 
main number. The quality of the structuring was checked via 
domain selective etching - more than 50 % of the desired 
domains are properly inverted. The radially poled part of 
the wafer was cut out and shaped (diamond-turned and pol- 
ished) into a WGR with the desired R/r ratio. The off-center 
parameter SR was estimated to be within 50 /im leading to 
SR/R w 0.03. 

Optical characterization of the fabricated WGR via line 
width and free-spectral-range measurements has shown that 
the intrinsic quality factor of the fundamental modes at A ~ 
1.04 /im is Q p « 4 x 10 7 and the large radius R ~ 1.58 mm. 
The last number is about 2.5 % larger compared to the de- 
sired one (1.54 mm); the difference is within the accuracy of 
our fabrication procedure. Correspondingly, phase matching 
is expected for j > 350 instead of 342. 

The setup for our optical experiments, depicted in Fig. 3, 
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FIG. 2. Tuning curves A S , ; (T) for A p = 1.04 /im, 2N = 686, 
R = 1.54 mm, r = 0.6 mm, and SR/R = 0, 0.003, and 0.03. The 
solid lines correspond to the strongest spectral peaks with \Sj\ > 
0.3, the dashed lines to 0.2 < \S 3 \ < 0.3, and the dotted lines to 
0.1 < \Sj\ < 0.2. The curves corresponding to \Sj\ < 0.1 are 
dropped. 
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FIG. 3. Schematic of the WGR-OPO experimental setup. The fre- 
quencies of the interacting waves are i/ p (pump), u s (signal), and Pi 
(idler). 



fled with the sum frequencies v v + v s and v v + v x , where u St j 
are the signal and idler frequencies, such that v s + v\ = u p . 
The accuracy of such measurements of ^j. s is within the line 
width. In some cases, not only a single process, but also two 
parametric processes, corresponding to adjacent values of j, 
are observed simultaneously. While no direct measurements 
of the signal and idler waves were possible with our spectrom- 
eter, the data obtained give a complete picture of the actual 
nonlinear processes. 

Changing the temperature T between 20 and 62 °C, we 
were able to tune the signal and idler wavelengths from 1.78 
to 2.5 ^m, including the point of degeneracy A s = Aj, in steps 
of (1 — 2) nm. The corresponding experimental data are pre- 
sented by the dots in Fig. 5. Most of the experimental dots 



comprises the whispering gallery resonator mounted on a 
heatable post. The latter one is used to change the resonator 
temperature up to 62 °C. An extraordinarily polarized pump 
beam at 1.04 /xm from an external-cavity diode laser is fo- 
cused into a rutile prism. In the focus, placed at the prism 
base, the pump beam couples into the WGR via frustrated to- 
tal internal reflection. In order to keep the setup as simple 
and as stable as possible, even under temperature expansion 
of the setup elements, the coupling prism contacts the res- 
onator rim. The residual pump wave and the generated waves 
coupled out of the resonator are guided to a spectrometer cov- 
ering the wavelength range from 0.2 to 1.1 fim. The spectra 
are collected for varying temperatures and pump powers. 

As soon as the pump wave is coupled into the WGR, two 
narrow peaks corresponding to the pump frequency v v and 
to the second harmonic 2v v become clearly visible in the 
spectrum. At pump powers P exceeding a threshold value 
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FIG. 4. Spectrum of the visible output of the WGR at a pump power 
of 10 mW for two different temperature values showing sum fre- 
quency peaks at v v + v s and v v + The corresponding parametric 
processes are phase matched at j = 351 and 352. 

P t h — 6 mW, parametric oscillation has been detected. In 
addition to the former spectral features, new peaks arise well 
above the noise level, see Fig. 4. They can be reliably identi- 
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FIG. 5. Signal and idler wavelengths (open and filled dots) evaluated 
from the frequency spectra as functions of the temperature T. The 
solid lines are the tuning curves calculated numerically for j = 350, 
351, and 352 from Eq. (2). 

nicely follow three tuning curves (solid lines) calculated for 
the neighboring spectral peaks with j = 350, 351, and 352. In 
accordance with theory, the measured temperature difference 
between two adjacent points of degeneracy is about 22 °C. 
Remarkably, filling of different tuning curves with the dots is 
non-uniform. For T < 30°C, the curve with j = 352 is prac- 
tically unfilled. The branches with j > 352 remain inactive in 
the whole temperature range. 

Several issues, closely related to the above results, are wor- 
thy of discussion: 

In some special cases, phase matching for \^ nonlinear 
processes can be realized in WGR's even without domain en- 
gineering, i.e. in the single-crystal case. With lithium niobate, 
it is possible for the processes involving modes of different 
polarizations fl3l EoTl . However, the tunability is restricted 
and the relevant nonlinear coefficient dia is about five times 
smaller than the largest coefficient ^33 B25I1 . 

Being focused on the tunability issue, we did not try to min- 
imize the pump threshold P t h- While our experimental value 
Pth ~ 6 mW is already low, it can still be decreased consid- 
erably. What are the prospects for decreasing the threshold 
power? In order to estimate them, we write down the scal- 
ing relation P t h oc ViV s /Q*Q*Q* d%, where dj 
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is the effective nonlinear coefficient and Q* s j are the loaded 
quality factors incorporating the coupling losses [26]. This 
relation contains the main parameters which can be varied in 
the experiment. Because of the imperfect domain engineer- 
ing, see above, the actual reduction factors \Sj\ « 0.1 were 
(3 — 4) times smaller than they might be, which gives roughly 
one order of magnitude increase in P t h- About two orders 
of magnitude increase stems from the overall coupling losses 
for p, s, i-waves - bringing the coupling prism in contact with 
the WGR rim noticeably decreases Q* s ; . Thus, the threshold 
power can be decreased by about three orders of magnitude 
without sacrificing the tuning properties. 

It is known that the phase-matching acceptance bandwidth 
increases nearby the point of degeneracy, A s = Aj, if the s, i- 
modes are of the same polarization [27]. This is why the ex- 
perimental dots in Fig. 5 are widely spread around the tuning 
curve for j = 351 at T « 45 °C. If a narrow acceptance 
bandwidth is required, one can employ an OPO scheme with 
different s, i-polarizations. The domain number 2N must be 
different in this case. About 1480 domains are needed, instead 
of former 2N rj 700, to realize this scheme for WGR radii 
similar to the ones in our experiment. The structure of the 
tuning curves is expected to be different as well. The above 
mentioned indicates a high degree of flexibility of the non- 
linear schemes based on structured WGR's. By changing the 
domain number one can proceed from broad-band to narrow- 
band gain and from identically polarized to cross-polarized 
generated waves. 

In conclusion, we have shown that radial poling of whis- 
pering-gallery resonators made of lithium niobate allows to 
combine a high tunability of nonlinear-optical processes, such 
as optical parametric oscillation, with common low-power 
continuous-wave light sources. The tuning characteristics dif- 
fer drastically from those known for bulk nonlinear schemes 
because of the discrete geometry-dependent mode structure. 
The potential of quasi-phase matching for shaping nonlinear 
processes is thus strongly extended. 

We thank the DFG and the Deutsche Telekom AG for fi- 
nancial support. 
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